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ABSTRACT 

We report on the results of a deep XMM-Newton observation of RBS 1774, 
the most recent dim isolated neutron star candidate found in the ROSAT archive 
data. Spectral and timing analysis of the high-quality PN and MOS data confirm 
the association of this source with an isolated neutron star. The spectrum is 
thermal and blackbody-like, and there is evidence at a significance level > 4cr 
that the source is an X-ray pulsar, with spin period of 9.437 s. Spectral fitting 
reveils the presence of an absorption feature at ~ 0.7 keV, but at this level data do 
not have enough resolution to allow us to discriminate between an absorption line 
or an edge. We compare the newly measured properties of RBS 1774 with those 
of other known dim isolated neutron stars, and discuss possible interpretations 
for the absorption feature. 

Subject headings: stars: individual (1RXS J214303.7+065419/RBS 1774) - 
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1. Introduction 


Over the last decade ROSAT observations have led to the discovery of seven very soft X- 
ray sources with quite particular characteristics. The extreme values of the X ray-to-optical 
flux ratio (> 10 4 ), the low hydrogen column densities (nn ~ 10 20 cm -2 ) and the detection 
of pulsations in the range ~ 5 — 10 s in four objects (five with the present one) strongly 
favor the possibility that these sources are close-by, X-ray emitting, dim, isolated neutron 
stars (XDINSs; see e.g. Treves et al. 2000 and Haberl 2004 for reviews and Table 1 for a 
summary). 

XDINSs stand apart with respect to other known classes of bona fide isolated neutron 
stars detected at X-ray energies. All of them are radio-silent, they are not associated with 
supernova remnants and exhibit a soft X-ray spectrum without evidence of a power-law hard 
tail. The origin of their X-ray emission has not been fully clarified as yet. The recent mea¬ 
surement of relatively large proper motions in three sources, implying transverse velocities 
> 150 kms -1 (Kaplan, van Kerkwijk & Anderson 2002; Motch, Zavlin & Haberl 2003; Motch 
et al. 2005), makes it unlikely that these are old (> 10' yr) neutron stars accreting from the 
interstellar medium (see e.g. Treves et al. 2000). Most probably the X-rays arise from the 
cooling of younger objects with an inferred age of ps 10 5 -10 6 yr. 

XDINSs play a key role in compact objects astrophysics being the only sources in which 
we can have a clean view of the compact star surface, without contamination from magneto- 
spheric emission or emission from a binary companion or a supernova remnant. They appear 
to be truly “isolated 1 ' neutron stars, and only a small number of them has been detected so 
far. Detailed multi wavelength studies of the largest possible sample of XDINS candidates 
are therefore fundamental for tracking the evolutionary history of galactic neutron stars, and 
for shedding light on their thermal and magnetic surface properties. 

Spectral analysis carried out so far has convincingly shown that the broadband X- 
ray emission from XDINSs is well represented by a blaekbody continuum, challenging the 
predictions of conventional atmospheric models. The absence of spectral lines appears to 
be well established in the brightest XDINS, RX J1856.5-3754 (Drake et al. 2002; Bur- 
witz et al. 2003), while absorption features have been increasingly reported in the spec¬ 
tra of other sources, including RX J1308.6+2127, RX J0720.4-3125, RX J1605.3+3249, 
RX J0806.4-4123, RX J0420.0-5022 (Haberl et al. 2003; Haberl et al. 2004b; Van Kerk- 
wjik et al. 2004; Haberl et al. 2004a). The absorption features are quite broad, in contrast 
with the atomic spectral lines predicted by atmospheric models, and they appear at energies 
~ 200 — 500 eV. Their nature is still uncertain but the intriguing identification with a proton 
cyclotron resonance has been suggested. Should this interpretation prove correct, it makes 
XDINSs highly magnetized objects with fields of the order of, or exceeding, a few xlO 13 G. 
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In this respect, the case of RX J0720.4-3125 is particularly interesting: the magnetic held 
strength derived by line energy, assuming it is a proton cyclotron feature, is in good agree¬ 
ment with that implied by the spindown measure (? ~ 3 - 6 x 10 _14 ss _1 ; Zaire et al. 2002; 
Kaplan et al. 2002; Cropper et al. 2004). 

We are still at a stage where every newly discovered XDINS holds important infor¬ 
mation for understanding the properties of the whole class and may also show some pecu¬ 
liarities which import unique information. RBS 1774 (1RXS J214303.7+065419) has been 
the most recent XDINSs to be found (Zampieri et al. 2001). The source is listed in the 
ROSAT Bright Source catalog, and lies about 48' off-axis in a PSPC pointing of the BL Lac 
MSS 2143.4+0704. The ROSAT observation had a limited statistics (~ 500 net counts), 
but it had been sufficient to reveal a very soft X-ray spectrum, to which an absorbed black- 
body provides an acceptable fit (kT ~ 92eV, rifj ~ 4.6 x 10 2O cm -2 ). No modulation with 
amplitude > 30% has been found in the PSPC data. Optical follow-up observations have 
shown no plausible optical counterpart in the X-ray error box down to R ~ 23, implying an 
X-to-optical flux ratio > 10 3 . All these properties are very similar to those of the already 
established XDINSs and make RBS 1774 worth of further investigation. 

In this paper we report results from a recent XMM-Newton observation of RBS 1774 and 
compare the newly measured properties of this source with those of other known members 
of this class. 


2. X-ray Observations 

RBS 1774 was observed with XMM-Newton on 2004 May 31 for 30 ks (EPIC PN). All 
the three EPIC detectors were configured in small window mode. We processed the data 
using the XMM-Newton SASv6.0.0. We selected periods of low background, with the Good 
Time Interval (GTI) hie leaving an effective exposure of 23 ks (EPIC PN). The same GTI 
was used for the spectral and timing analysis. 

The source is clearly visible in EPIC PN, MOS1 and MOS2 with count rates of 2.054 + 
0.011, 0.346 + 0.004 and 0.328 + 0.004 counts s^ 1 respectively (0.2-2 keV band). In the EPIC 
image, the source is slightly off axis since the pointing was requested based on the revised 
position published by Zampieri et al. (2001). The position obtained with XMM-Newton 
(J2000 coordinates) is a = 2lM3 m 03.3 s , 6 = +06 o 54'17", with a 90% uncertainty radius of 
3". Within 3 cr, previous estimates based on ROSAT data (see again Zampieri et al. 2001) are 
consistent with the XMM position. In particular, the position from the ROSAT ASS Bright 
Sources Catalogue (a = 21 /l 43 m 03.7 s , 5 = +06°54'19.5", 90% uncertainty radius of 18"; 
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labelled with RASS in Zampieri et al. 2001) is within ~ 7" from the present XMM position. 
Figure 1 shows the ROSAT and XMM-Newton error boxes overlayed on the optical image 
(discussed below). We note the small difference in the position of the RASS error box with 
respect to that reported by Zampieri et al. (2001). 


2.1. Spectral analysis 

For the EPIC PN, we extracted single and double events (pattern < 4, FLAG=0) 
within a circle of 25" radius for source and within two circular regions (to avoid out-of- 
tirnc events) with radius 50" for background. For the EPIC MOS1 and 2, we extracted up 
to quadruple events (pattern < 12, FLAG=0) with a 25" and 32" radius for the source 
respectively with two rectangular background regions to exploit the limited window size more 
effectively. In order to account for the slightly off-axis position, we generated the appropriate 
off-axis response and ancillary response matrices. Data were grouped, ensuring a minimum 
of 30 counts per energy bin. 

We performed a spectral analysis by fitting data from the three EPIC detectors simul¬ 
taneously. We find that the broadband spectrum can be represented by a single absorbed 
blackbody, with kT = 0.1014 keV and nn = 3.65 x 10 20 cm -2 (see Figure 2, Table 2). Both 
parameters are consistent with those measured in the past using ROSAT data, although bet¬ 
ter constrained. The unabsorbed EPIC-PN flux in the 0.1 — 2.4 keV band is 6.05 x ICG 12 erg 
cm -2 s“\ slightly less than that reported by Zampieri et al. for ROSAT PSPC (8.7 x 10” 12 
erg cm -2 s -1 (the difference is likely to be due to the large off-axis angle of RBS 1774 in 
the ROSAT observation). Assuming a typical luminosity of 5 x 10 31 erg/s gives a distance 
~ 280 pc. The value of the column density is similar to, but slightly lower than, the total 
Galactic absorption in the source direction (tih = 5 x 10 20 cm -2 , Dickey & Lockman 1990). 
The resulting reduced y 2 = 1.36 is not fully satisfactory and an inspection of the residuals 
suggests that the largest discrepancies between model and data are above ~ 0.6 keV. 

As discussed by many authors (see e.g. Zavlin & Pavlov 2001, Treves et al. 2000 for 
reviews), realistic model spectra of the cooling atmosphere surrounding an unmagnetized 
(B < 10 9 G) neutron star exhibit a distinctive hardening with respect to a blackbody 
and deviate significantly from a Planckian shape. When accounting for a magnetic field of 
moderate strength (B ~ 10 12 — 10 13 G), the hard tail present in nonmagnetic models with 
comparable luminosity is partially suppressed and the X-ray spectrum, although still harder 
than a blackbody at the neutron star effective temperature, T e //, is more Planckian in shape. 

Models of fully ionized, pure H neutron star cooling atmospheres, as computed by 
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Pavlov et al. (1992) and Zavlin et al. (1996), are currently implemented in XSPEC (nsa 
model) for three different values of the magnetic field: B — 0 (unmagnetized models) and 
B = 10 12 ,10 13 G. In the attempt to improve the spectral fit we tested all the three sets 
of models, but in all cases the resulting y 2 was worse than that obtained using a simple 
blackbody model (see Table 3). The star mass and radius have been fixed at M = 1.4M 0 
and R = 10 km, but the effects of gravity on the emergent spectra are too small to produce 
any appreciable difference in the fits. 

Similarly, no improvement is obtained adding a power-law component or a second black- 
body. However, we find that adding a gaussian line in absorption at ~ 0.7 keV produces 
a statistically significant (Ay 2 = 56; 7a) improvement 1 , leading to a reduced y 2 = 1.20. 
The best fitting energy of the line is Eu ne = 0.754 keV; the line width and depth are 
o~Hne = 0.027 keV and ru ne = 4.8, respectively. The latter two parameters are mutually de¬ 
generate and not well constrained, because, at least for narrow or moderately narrow lines, 
the energy resolution is insufficient to strongly distinguish between width and depth. A 
slightly better (reduced y 2 = 1.17) but similar fit is obtained by using an absorption edge 
instead of a gaussian line (see Figure 3): in the edge model there is one parameter less, there¬ 
fore the depth is well constrained. At this level the data do not have enough resolution to 
allow us to discriminate between the two kinds of spectral features, in terms of their physical 
meaning. The best fitting energy of the edge is E e d ge = 0.694 keV and the optical depth is 
Tedge = 0.25. In both cases, the blackbody temperature is essentially unchanged with respect 
to the single blackbody fit. The best fit parameters are reported in Table 2. The flux in 
the line/edge has been computed by fitting data with the models above, then switching the 
optical depth to zero and calculating the difference in flux, yielding 8.6 x 10~ 14 erg cm _2 s _1 
(edge) and 3.3 x 10” 14 erg cm _2 s _1 (gaussian line). 


2.2. Timing analysis 

The high throughput of XMM-Newton provides a much more stringent check than the 
ROSAT PSPC on whether RBS 1774 is pulsating. For the timing analysis, we extracted the 
counts within a 30" radius aperture in both PN and MOS1 and MOS2 but with the aperture 
truncated by a chord where it approached the edge of the CCD window in each case: this 
maximised the counts while avoiding any edge effects. A background area was selected using 


1 The F-test has been widely used to test the significance of spectral lines, although this is strictly in¬ 
appropriate (see Protassov et al. 2002). Only for completeness, we report that we have checked the F-test 
statistics value and probability, obtaining 12.1 and 1.8 x 10 7 , respectively. 
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as much as possible of the small window mode image area while avoiding the edge of the 
window and the out-of-time events. This resulted in event hies containing 38485, 9738 and 
9500 events available for data analysis for PN, MOS1 and MOS2 respectively. 

We analysed PN, MOS1 and MOS2 data both individually and combined using the 
Maximum Likelihood Periodogram technique described in Zaire et al. (2002) and Cropper 
et al. (2004). We searched periods from 10000 s to 30 ms, ensuring that in each case the 
period grid was 2.5 times better sampled than the Nyquist frequency. Below 0.3 s the time 
resolution even in small window mode of MOS prevents them being used, so only the PN 
data were used for periods shorter than this. 

The search revealed a significant period at 9.437 s. This is illustrated for the combined 
MOS and PN dataset in Figure 4, where the dotted lines show confidence intervals of 68% 
and 90%. The Ay 2 between the maximum and the noise floor is 40 in the combined dataset, 
while for the PN alone Ay 2 = 36. The peak is evident in the combined MOS1 and MOS2 
data, with Ay 2 = 9.5. We have checked whether the reduced significance of the MOS1+2 
peak is consistent with that in the PN data, given its greater sensitivity, by selecting (at 
random) event data from the PN event list to generate a reduced dataset with the same 
number of events as MOS1+2. In this reduced PN dataset the noise levels are similar to 
the MOS1+2 dataset, and the 9.437 s peak is no longer the most significant in the 9-10 sec 
period range. This indicates that the lower significance of the peak in the MOS1+2 dataset 
is consistent with its reduced statistics with respect to the PN dataset and that no additional 
bias is introduced because of event selection. 

Ay 2 = 40 corresponds to a > 6a result (a = \J Ay 2 for 1 degree of freedom). However, 
we need to consider the random probability of a peak occurring given the large number of 
periods searched. In the interval between 1 and 1000 seconds, there are 59638 independent 
periods, and the probability of a 6er peak will occur once in 1.5 x 10 5 cases, corresponding to 
a Ay 2 = 18.7 or > 4 a result. We therefore consider the peak to be significant, although with 
the caveat that in Haberl, Pietsch & Motch (1999) a 4a detection was claimed for a 22.69 s 
periodicity in RX J0420.0-5022, which has since been shown to have a 3.45 s periodicity. 

We have consequently folded the data on the 9.437 s period and binned them in phase. 
This is shown in Figure 5 for PN separately, for the combined MOS data, and then for all 
three instruments together. The MOS data taken separately are noisy, with a variation which 
would not be considered significant in its own right, but which nevertheless shows an ampli¬ 
tude and phasing which is consistent with that from the PN. This is evident (equivalently) 
from our periodogram analysis where the Ay 2 increases from 36 to 40 when the MOS1+2 
data are added to the PN data, and adds weight to our conclusion that this pulsation is real. 
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The pulsation is approximately sinusoidal, and the amplitude of the variation is remark¬ 
ably small: fitting a sinusoid yields a semiamplitude of 0.036 ± 0.006 (0.072 peak-to-peak). 
In the attempt to confirm the periodicity, we re-analysed the archival ROSAT data taken in 
May 1991 (Zampieri et al. 2001). This observation is divided into two segments of duration 
~ 1800 s each; the interval between the two segments is ~ 4000 s. No statistically significant 
peak is present in the power density spectrum of the two separate segments, due to the poor 
statistics. No modulation with amplitude > 30% has been found in the PSPC data, and this 
limit is consistent with the XMM-Newton detection. 


3. Possible nature of the spectral feature 

The quality of the blaekbody fit to the phase-averaged EPIC data of RBS 1774 is not 
completely satisfactory. Adding an absorption edge or a gaussian line in absorption signifi¬ 
cantly improves the spectral fit, suggesting the presence of an absorption feature similar to 
those observed in RX J1308.6+2127, RX J0720.4-3125, RX J1605.3+3249, RX J0806.4-4123, 
RX J0420.0-5022 (Haberl et al. 2003; Haberl et al. 2004b; Van Kerkwjik et al. 2004; Haberl 
et al. 2004a). However, the absorption features detected in these other sources appear at 
energies below 500 eV, while in the case of RBS 1774 it is at ~ 700 eV (for the purpose 
of the discussion we use the energy of the feature obtained by the blaekbody plus edge fit, 
which is slightly better constrained). This makes the detection in RBS 1774 less sensitive to 
the uncertainties in the EPIC-PN spectral calibration, since residual systematic calibration 
effects are mainly present below ~ 0.5 keV. We note that a feature at this energy has been 
detected in the spectrum of the neutron star IE 1207.4-5209 (Sanwal et al. 2002; Mereghetti 
et al. 2002a; Bignami et al. 2004; Hailey & Mori 2002; Mori et al. 2004; Zavlin et al. 2004). 
The comparison between the properties of the XDINSs and those of IE 1207.4-5209 must 
be done carefully, given the fact that the latter is probably younger: it has a shorter period, 
it is still associated with a supernova remant and it shows an additional harder spectral 
component in the X-ray spectrum. Also, IE 1207.4-5209 is unique since it exhibits at least 
two spectral lines, with harmonic spacing. However, as will be discussed further on, some of 
the models proposed for this source can be applied to RBS 1774 as well. 

Strong atomic absorption features are predicted by low-field atmospheric models with 
significant metal abundances (either solar, or pure iron, see Rajagopal & Romani 1996; 
Zavlin et al. 1996). The predictions for strongly magnetized metal-enriched atmospheres 
are more difficult to assess, due to the lack of detailed calculations. However, the strongly 
magnetized pure Fe spectra presented by Rajagopal, Romani & Miller (1997) also show a 
significant number of absorption lines. The main difficulties faced by these models is that 



they exhibit a variety of rather narrow lines, in contrast with the single feature observed in 
RBS 1774. 


Other possible interpretations are discussed below and can be divided in two groups, 
requiring either a superstrong (B ~ 10 14 G) or a moderate (B ~ 10 10 — 10 12 G) magnetic 
held. 


3.1. Super-strong magnetic fields 

A plausible hypothesis is that the feature in the spectrum of RBS 1774 is associated with 
a proton cyclotron line. An interpretation in terms of a proton resonance has been recently 
put forward for the lines observed in other XDINSs, and, at least for RX J0720.4-3125, it is 
supported by the agreement between the magnetic held strength independently inferred from 
spin-down measurements and line energy (Cropper et al. 2004). In the case of RBS 1774, 
this would imply an ultra-strong magnetic held of B ~ (E edge / 0.63keV)(l + z) x 10 14 G ~ 
1.4 x 10 14 G, assuming a gravitational redshift factor at the star surface of (1 + z)^ 1 = 
y/1 - 2GM/c 2 R 0.8. Such a held strength is comparable (although a factor ~ 30% 
higher) with that estimated in a similar way for RX J1605.3+3249 and RX J0806.4-4123 
(~ 9 x 10 13 G, see Van Kerkwjik et al. 2004; Haberl et al. 2004a), and is well above the 
critical value, Bq ~ 4.4 x 10 13 G, at which QED effects become important. 

An alternative interpretation is that we are detecting an absorption feature due to 
a bound-bound or bound-free transition in the hydrogen atom. A strong magnetic held 
can signihcantly increase the atomic binding energies and radiative transition rates, mov¬ 
ing the H ionization edge into the soft X-ray band. First calculations of hydrogen at¬ 
mospheres with bound atoms have shown that broad spectral lines become prominent if 
the atmospheric temperature drops below ~ 10 6 K. The strongest line corresponds to 
the transition between the ground state and the lowest excited state at an energy of « 
75[1 + 0.13ln(B/10 13 G)] + 63(B/10 13 G) eV (Zavlin & Pavlov 2001). If applied to the 
~ 0.7 keV line detected in RBS 1774, this interpretation yields a magnetic held of B m 10 14 G, 
again in the ultra-magnetized range. More detailed atmospheric models which include par¬ 
tially ionized hydrogen have been computed by Ho et al. (2003), accounting for magnetic 
held effects and using the latest available equation of state. They have shown that, for 
B = 5 x 10 14 G and T = 5 x 10 6 K, two broad absorption features, due to bound-free 
transitions to different continuum states, are present at ~ 0.76 and 4 keV, respectively. The 
latter is much stronger, but it would not be visible in the spectrum of RBS 1774 which 
has negligible emission above ~ 2 keV. At somewhat lower helds (B ~ 10 14 G) Ho & Lai 
(2004) found that the dominant atomic features in the 0.2-2 keV band are those due to 
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the low-energy bound-free transition, which is barely visible at 0.54 keV, and the bound- 
bound transition discussed above at 0.74 keV. These authors argue that at super-strong 
held strengths (B > 10 14 G) vacuum polarization may efficiently suppress these lines, in 
which case the above conclusions can be (at least partially) invalidated. However, current 
treatment of this QED effect is too crude to make a definite statement. 

At present, accurate atmospheric models for chemical compositions other than hydrogen 
are unavailable for B > 10 12 G. In the case of IE 1207.4-5209, attempts to interpret the 
~ 0.7 keV feature as a transition in ions other than H have been made, simply based on the 
latest theoretical results on the characteristics of bound-bound and bound-free transitions in 
He-like and H-likc ions in strong fields. In this case, a ~ 0.7 keV feature may be associated 
with a transition between the ground level (m = 0,^ = 0, where m and v are the magnetic 
and longitudinal quantum numbers) and the lowest excited tightly-bound level (m = 0, v — 
1) of He + . The required magnetic held is ~ 1.4 x 10 14 G, for (1 + z)~ l ~ 0.81 (Sanwal 
et al. 2002; Zavlin et al. 2004). Taken overall, these results show that if RBS 1774 has a 
magnetar-like held, B « 10 14 G, the feature at ~ 0.7 keV may be due to a proton cyclotron 
resonance, or bound-bound transitions in H or H-like He, or a blend of these. 


3.2. Moderate magnetic fields 

An obvious scenario involving a moderate magnetic field is that in which the observed 
feature is associated with electron (instead of proton) cyclotron resonance. Would this be 
the case, the surface magnetic field is necessarily much lower, B ~ (E edge /11.6 keV)(l + z) x 
10 12 G ~ 7.5 x 10 10 G for (1 + z)- 1 ~ 0.8. 

A further possibility is that of an atomic transition out of the (m = 0, v = 0) ground 
state to the (m = 0, v = 1) excited state or to the continuum in a rnid-Z H-like ion, such 
as C VI, N VII or O VIII. Such an interpretation has been put forward by Hailey & Mori 
(2002) again in connection with the lines in IE 1207.4-5209 and yields B ~ 10 11 G (we refer 
in particular to their solutions B; see also their discussion about the possibility of having 
elements other than H in the neutron star atmosphere). Both these models are conceivable, 
although the required field strength is well below that estimated for RX J0720.4-3125 and 
well below the average fields measured in ordinary radio pulsars (only ~ 20% of the pulsar 
population has B < 6 x 10 11 G). 
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4. Comparison with other thermally emitting Neutron Stars and open issues 

In this paper we have reported on the results of a recent XMM-Newton observation of 
RBS 1774. Spectral and timing analysis of the high-quality PN and MOS data confirm the 
association of this source with an isolated neutron star (as originally reported by Zampieri 
et al. 2001). The spectrum is thermal and there is evidence at a significance level > 4cr that 
the source is an X-ray pulsar, with spin period of 9.437 s. The EPIC-PN X-ray light curve 
is single peaked and approximately sinusoidal, in agreement with what is seen in most of 
the other XDINSs (see Cropper et al. 2001; Haberl et al. 2004a; RX J1308.6+2127 is the 
only XDINSs with double peaked X-ray pulse). However, both the relatively low level of 
significance of the periodicity and the small amplitude variation in the folded data prevented 
us from performing pulsed phase spectroscopy. 

The deepest available optical pointing of the field is that obtained with the New Tech¬ 
nology Telescope (NTT) at la Silla on 2001 May 27. The faintest source in this image is 
object D, as listed in Table 2 by Zampieri et al. (2001), with R = 22.77. This, combined 
with the unabsorbed XMM-Newton X-ray flux (f x = 5 x 10“ 12 erg cm _2 s _1 in the 0.2-2 keV 
band for Tin = 3.6 x 10~ 20 cm -2 ) allows a lower limit to be placed on the X-to-optical flux 
ratio of fx/fn 1700. When using the absorbed flux f x = 2.8 x 10 12 erg cm 2 s 1 and 
the same X-ray energy band, we get fx/fn ~ 1000. 

Taken overall, the characteristics of RBS 1774 are remarkably similar to those of an 
other XDINS recently studied in detail, RX J0806.4-4123 (see Haberl et al. 2004a). Both 
sources have a similar EPIC-PN count rate of ~ 1.7 — 2.0 counts s -1 and a similar period 
(~ 10 s). 2 Also, they have a blaekbody temperature of ~ 100 eV, among the highest observed 
from the XDINs, and exhibit the smallest pulsed fractions of the group (semiamplitude of 
4% and 6% respectively) which makes the detection of significant spectral variations with 
phase difficult. Intriguingly, and if related, the two latter characteristics may be suggestive 
of a nearly pole-on viewing geometry, with the spin and magnetic axes fairly closely aligned. 
This appears also broadly consistent with the fact that the energy of the spectral feature is 
the highest in RBS 1774, RX J0806.4-4123 and RX J1605.3+3249. If the feature energy is 
sensitive to the field strength, then, both the temperature and the average magnetic field 
are higher in the polar regions. 

As discussed in § 3, many of the viable interpretations for the feature detected at 
~ 0.7 keV require a magnetar-like field strength, B ~ 10 14 G. For a period of 9.437 s, this 
implies a dipole spin-down rate (at present) of P ~ 11 x 10~ 13 s/s, and a spin-down age of 


2 RX J1308.6+2127 also has a similar count rate, but a much higher pulsed fraction. 
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1.4 x 10 5 yrs (assuming B constant). In order to investigate if such a high magnetic held 
strength is compatible with the cooling age of the source, we considered the evolutionary 
tracks in the B-P diagram by repeating the simple computation discussed in Cropper et al. 
(2004). We computed the source age, Td, and the value of the magnetic held at the birth 
of the neutron star, Bo, obtaining Bo ps 1.6 x 10 15 G and Td ~ 10 4 yrs for Hall cascade; 
Bo ss (1.1 or 1.5) x 10 14 G and Td ~ (1.4 or 1.1) x 10 5 yrs for ambipolar diffusion in the 
solenoidal or irrotational mode, respectively. The situation is therefore very similar to that of 
RX J0720.4-3125 (Cropper et al. 2004): ambipolar diffusion predicts a magnetic field which 
is relatively constant over the source lifetime, while Hall cascade gives a scenario compatible 
with a superstrong field at the star’s birth ( B 0 ~ 10 15 G). The corresponding star’s ages are 
compatible with the cooling age of a neutron star with mass ps 1.35 M 0 or slightly larger, 
1.4 — 1.5 M 0 (see Fig. 2 in Yakovlev et al. 2004). 

Should the strength of the magnetic field be confirmed by future observations, this will 
make RBS 1774 the XDINS with the highest field discovered so far. RX J1605.3+3249 and 
RX J0806.4-4123 may possess magnetic fields of comparable strength (B ~ 9 x 10 13 G Van 
Kerkwjik et al. 2004; Haberl et al. 2004a), and this raises the question of whether or not a 
highly magnetized neutron star may remain for long periods without any sign of some kind 
of transient activity. Very recently, a few radio pulsars with B > Bq and long period (~ few 
seconds) have been discovered (see e.g. McLaughlin et al. 2003; Morris et al. 2002; Camilo 
et al. 2000). On the other hand, non-thcrmal magnetospheric emission from XDINSs is so 
far undetected, although this may simply reflect the fact that it falls below the detection 
threshold of present instruments. In fact, if XDINSs are similar to the standard radio¬ 
pulsars the relation between rotational energy and non-thermal X-ray luminosity would be 
L\ ~ 10~ 3 if. Taking the above value for the magnetic field gives E ps 10 31 erg/s, for 
RBS 1774, implying L x ps 10 28 erg/s. 

XDINSs were unanimously believed to be steady sources, as indicated by several years 
of observations for the brightest of them. Only recently, and quite surprisingly, XMM- 
Newton observations revealed a substantial change in the spectral shape and pulse profile of 
RX J0720.4-3125 over a timescale of ~ 2 yr (de Vries et al. 2004; Vink et al. 2004). Possible 
variations in the pulse profile of RX J0420.0-5022 over a similar timescale (~ 0.5 yr) have also 
been reported (Haberl et al. 2004a), although at a much lower significance level. Even more 
recently, timing analysis of Chandra and XMM-Newton data of IE 1207.4-5209 revealed non¬ 
monotonic variations in the spin period, possibly related to glitches occurring on a timescale 
of ~ 1-2 yr (Zavlin et al. 2004). Also the monitoring of anomalous X-ray pulsars (AXPs), 
previously believed to be steady emitters, has now revealed X-ray variability: both short 
energetic bursts and long-lasting variations in the X-ray emission of some of these sources 
have been reported (see Gavriil, Kaspi & Woods 2003, Woods et al. 2004, and Mereghetti 
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et al. 2004 for the cases of IE 2259+58 and IE 1048-59). Quite interestingly, the long¬ 
term changes observed in RX J0720.4-3125 bear some resemblance to those reported in IE 
2259+58 and IE 1048-59, with the difference that in RX J0720.4-3125 the spectral and 
lightcurve evolution is not accompanied by any flux variation. 3 

Further timing studies of XDINSs will be of key importance in securing new spin-down 
measurements: this will provide a second independent measurement of the held strength and 
also will bring further insights on the long-term behavior of XDINSs. Should the presence of 
ultra-strong fields and the variability be confirmed, the putative relation between XDINSs 
and AXPs, put forward on the basis of the similarity in spin periods (see e.g. Mereghetti et al. 
2002b, Haberl 2004), becomes much firmer. In turn, this relation may provide the “missing 
link” among XDINSs and SGRs. If indeed XDINSs are the descendants of SGRs/AXPs, 
as their lower temperature may suggest, then explaining why evolution produces “quiet” 
magnetars is still a challenge. 


We thank G.G. Pavlov for a discussion about He + transitions in the presence of strong 
magnetic fields, in connection with the model presented by Sanwal et al. (2002). RT, LZ 
and AT acknowledge financial support from the Italian Ministry for Education, University 
and Research through grant PRIN-2002-027245 and PRIN-2004-023189. JJD acknowledges 
support from NASA contract NAS8-39073 to the Chandra X-ray Center. We thank an 
anonymous referee for several useful comments on the manuscript. 
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Table 1. Summary of XDINSs properties 


Source 

UT°° 

K1 bb 

eV 

n H 

10 2 ° cm -2 

Pline 

eV 

P Semi-Ampl. a 

s 

Ref. b 

RX J1856-3754 

56.7 

0.18 

no 

no 

no 

1 

RX J0720-3125 

85.2 

1.38 

270 

8.39 

n% 

2 

RX J1605.3+3249 

94.1 

0.68 

493 

no 

no 

3 

RX J1308.6+2127 

85.8 

4.10 

290 

10.31 

18% 

4 

RX J0420-5022 

44.9 

1.02 

329 

3.45 

13% 

5 

RX J0806-4123 

95.6 

0.41 

460 

11.37 

6% 

5 

RBS 1774 

101.4 

3.65 

700 

9.44 

4% 

6 


a Semi amplitude of the folded light curve. 

b l) Burwitz et al. (2003); 2) Haberl et al. (2004b); 3) Van Kerkwjik et al. 
(2004); 4) Haberl et al. (2003); 5) Haberl et al. (2004a); 6) this paper. 


Table 2. Model Fit Parameters: blackbody (bb) model. a 


Model 

n H 

K1 bb 

rp 

edge/line 

'Tedgei 1~line 

&line 

fx h 

y 2 /d.o.f. 


10 2U cm -2 

eV 

eV 


eV 

erg cnr“ 2 s _1 


bb 

Q £C+0.16 
O.OO—o.i3 

101.4±g;| 




5.16 x 10” 12 

1.36 

bb+abs. edge 

q 6 n+ 0 - 21 

o-OU-0 16 

104.0+D;? 

6941^ 

0 2 C > +0 - 03 
u - zo -0.03 


5.07 x 10" 12 

1.17 

bb+ gauss, line 

o 74+°- 14 

o- '^-0.10 

io2.it8:i 

754lg 

4 '^-0.5 

27lg 5 

5.20 x 10" 12 

1.20 


a All fits have been obtained by fitting simultaneously data from EPIC PN, MOS1 and MOS2; the 
absorption model is TBabs in XSPEC. The reported errors are the 68% (ler) confidence range. The 
number of degrees of freedom (d.o.f.) is: 286 (bb), 284 (bb+abs. edge) and 283 (bb+gauss. line), in the 
three cases respectively. 

b Unabsorbed flux measured with EPIC-PN in the (0.2 — 2) keV band. 











Table 3. Model Fit Parameters: atmospheric models 3. 


B 

n H 

kT°° 

X 2 /d.o.f. 

10 12 G 

10 20 cm -2 

eV 


0 

8%±a2 

29-8±8;| 

1.94 

1 

o q+0.2 

47.3lgj 

1.99 

10 

8.6183 

50.0183 

2.19 


a All fits have been obtained by fitting si¬ 
multaneously data from EPIC PN, MOS1 
and MOS2; the absorption model is TBabs 
in XSPEC. The reported errors are the 
68% (lcr) confidence range, and the num¬ 
ber of degrees of freedom (d.o.f.) is 286 
in all the three cases. The star mass 
and radius have been fixed at 1.4 M 0 
and 10 km, respectively. The tempera¬ 
ture measured by a distant observer, T°°, 
and the effective temperature obtained by 
the model fit, T e //, are related by T°° = 
(1 + z)~ 1 T e f f , where (1 + z)- 1 = [1 - 
2.95(M/M 0 )/(i?/lOkm)] 1 / 2 is the gravita¬ 
tional redshift factor. 
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Fig. 1.— R-band NTT image of the field around RBS 1774 (limiting magnitude ~ 22.8, see 
text and Zampieri et al. 2001). The red circle (3” radius) shows the source position as derived 
from our recent XMM-Newton observation, while the blue circles show the position inferred 
from past data (see the electronic version for a full color image). Positions (see Zampieri et 
al. 2001 for all details) are taken from: White et al. (1994) Catalogue (WGA); ROSAT ASS 
Bright Sources Catalogue (RASS, Voges et al. 1996), ROSAT SRC Catalogue (ROSATSRC, 
Zimmermann 1994), RASS held analized with a Brera Multiscale Wavelet alghoritm (BMW, 
Lazzati et al. 1999). Error circles represent the 90% confidence intervals. 
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Fig. 2.— The top panel shows the count rate spectrum of RBS 1774 obtained with EPIC 
PN (black) and the two EPIC MOS (light green, red in the color version) detectors, with the 
best single blackbody combined model fit (parameters are given in Table 2). The bottom 
panel shows the data to model ratio. 
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Fig. 3.— Top panel: the count rate spectrum of RBS 1774 obtained with EPIC PN and the 
two EPIC MOS detectors. The best fit model is an absorbed blackbody with an absorption 
edge at ~ 0.7 keV (parameters are given in Table 2). The bottom panel shows the data to 
model ratio. 





































































